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Abstract 

We study the forced rupture of adhesive contacts between monomers that are not covalently hnked 
in a Rouse chain. When the appUed force (/) to the chain end is less than the critical force for 
rupture {fc), the reversible rupture process is coupled to the internal Rouse modes. If ///c>l 
the rupture is irreversible. In both limits, the non-exponential distribution of contact lifetimes, 
which depends sensitively on the location of the contact, follows the double-exponential (Gumbcl) 
distribution. When two contacts are well separated along the chain, the rate limiting step in the 
sequential rupture kinetics is the disruption of the contact that is in the chain interior. If the 
two contacts are close to each other, they cooperate to sustain the stress, which results in an 
"all-or-none" transition. 



PACS numbers: 



Nanomanipulation of biological molecules by force can be used to explore and control 
their energy landscapes at the single molecule level [Tj. The challenge is to solve the inverse 
problem, namely, to extract the unbinding pathways and structural features pertaining to 
the systems of interest from measurements such as the force-extension curves, lifetime dis- 
tribution of adhesive contacts, and unbinding force distributions. From these data, many 
characteristics of the energy landscape (roughness, barrier height, and movement of the 
transition states) have been obtained using theoretical models [3] . In these models, the 
force-induced unbinding transitions from the bound state (for example folded conformations 
of RNA) are presumed to occur irreverisbly along a one dimensional reaction coordinate, 
R, that is conjugate to the applied force, /. In this picture, / "tilts" the energy landscape, 
thus rendering the bound state unstable when / exceeds a critical value. On the other 
hand, detailed simulations of forced unfolding of proteins and RNA ^ show that collective 
topology-dependent events that describe the internal dynamics [3] are coupled to R. Thus, 
it is important to develop solvable models for force-induced disruption of "folded" structures 
starting from a well-defined Hamiltonian. 

In this letter, we illustrate the interplay between internal polymer dynamics and the 
rupture kinetics of intramolecular adhesive contacts between monomers. Such a model is a 
caricature of forced-unfolding of RNA hairpins [1] where the adhesive contact mimics the 
stability of the stem. Using an exactly solvable Rouse model, with one or two adhesive 
contacts between monomers that are not covalently linked, we address the following ques- 
tions: (i) For a single adhesive contact, how do the lifetime distribution functions vary with 
the backbone distance between monomers r and s that are in adhesive contact? (ii) To 
what extent is the binding of a ruptured contact influenced by / and the internal polymer 
motions? 

The model: Consider a Rouse chain with identical monomers, connected by a harmonic 
spring with the equilibrium distance a, with a single adhesive contact (r, s) between two non- 
covalently linked monomers r and s, l<r<s<N. The first monomer is fixed, and a constant 
force f=/x is applied to the A^*'^ monomer in the direction parallel to the end-to-end vector 
X. Letting Ilrn={Rm^ ^m} ^e the monomer position, {m=l,2, . . . , N), the equation 
of motion is, ^J:Rm=— VR^i!f({Rm})+Fm where VR„=9/(9Rm, ^ is the friction cefficient, 
and Fm is the random force with zero mean and {F^(t)F^{0))=2C^kBT6mnSai3S(t) (a, l3=x, y, 
z). We model the adhesive contact (r, s) using a harmonic potential with a spring constant. 
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K,=3kBT/b'^, where b is the equihbrium contact distance. The Hamiltonian with a single 
contact under tension is taken to be H=HQ+f-yi{Il]\r — Ri) where 

^0 = ^ E(R. - Rn-i)^ + ^e(5 - |R..|)R^, (1) 

n=2 

with Rrs=Rr.— Rs, and Q{x) is the Heaviside step function. The contact is ruptured when 
IRtsI^-B. The probabihty density function (pdf) of the (r, s)-contact hfetimes is obtained as 
Prs{t, f)=—dPrs{t, f)/dt, where Prs{t), the probabihty that the adhesive contact remains in- 
tact, IS Prsit, f)= 7rf3R^^jrf3R0^Q^(R^^^i.R0jp^(R0j • Here, Q/(R.„t; R^J, the conditional 
probability for the contact distance 'R,rs with the initial distribution P/(R[!g) of the contact 
distance, R^^, is given in terms of Ci(t) = (Rr.s(t)), and C2(t) = ([Rrs(t) — (Rrs(t))]^) by 

We define the A^x3-matrices for positions, R(t)=(Ri(t), R2(t), . . . , RAr(t))^ random forces, 
F(t) =(Fi(t), F2(t), . . . , F7v(t))^ and the applied force, f=(0, 0, . . . , fJv)^ where "f" denotes 
the adjoint 3xA^ matrix. The equation of motion in matrix form is J:R + MR = ^F + 

with the NxN matrix M=co'o(Ao+7Ai). The tridiagonal Rouse matrix is Aq, with 
«i,i=«Ar,Af=l' «m,m=2, and ali^,m+i=am_i,m=- 1 [5], and uJo=3kBT/2a'^. The interaction 
between monomers r and s are included in Ai, with al^=alg=l and 0^:^=0^^=— 1. The 
strength of the adhesive contact, relative to the harmonic bond is 7=(a/6)^. 

Integration of the equation of motion gives R(t)=G(t)R(0) + (l/0/o c?i'G(f)(F(t')+f ), 
where G(t) is the Greens function; G{uj) = [iujl + M]^^ in the Fourier domain. 
To compute Ci(t) and C2(t), we must invert G(ci;). We rewrite Ai=PP"l", 
where P is an A^xl vector with non-zero elements Pr=l and Ps=—1, and 
G(w)=Go(cj)-Go(cu)P[;^I + PtGo(cj)P]-^PtGo(cu) |], where the matrix elements 
Go{uU=[zuI+Mo];A=j^Ep and X,=nV/N' are the eigenmodes 

(p=o,i,...,iv) [g. 

The pdf of contact lifetimes for a single contact: Constraining the first residue of the 
chain is equivalent to excluding Ao=0, which defines the overall translation of the chain. The 
results are presented by including the first ten modes for A^=100, a=OAnm, the diffusion 
coefficient, D=10^nm'^/s, and kBT=4.1pNnm. We chose Pj(RQJ = ^^^^^^^, (5(|R°J - B/2), 
with the rupture distance set to B=2a. To assess the effect of the strength of the adhesive 
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FIG. 1: The survival probabilities Prs{t) for a single contact (r, s)=(l, 100) (a, d), (20,80) (b, 
e), and (40,60) (c, f), with /=lpN (a-c) and 40pN (d-f), and with 7=1 (solid lines) and 0.01 
(dashed lines). The lifetime distributions Prs{t) are shown in the insets. Notice the variations in 

time scales in the panels. The chain structures for the intact single contact (20,80), before and 
after the rupture are shown. 

contact on the contact rupture kinetics, we used 7=1 (strong contact) and 7=0.01 (weak 
contact). The critical equilibrium force required to rupture the contact is /c=(3A;bT_B/6^); 
fc—31pN for 7=1, and fc—0.31pN for 7=0.01. The values of / were set to f—lpN (weak 
force) and f^AQpN (strong force). At f^lpN, ///c<l for 7=1, and ///c~l for 7=0.01; 
at /=40p7V, ///e>l for 7=1 and ///c>l for 7=0.01. 

The survival probabilities Prs(^) (Fig- 1) of the adhesive contact at f=lpN and f=AOpN 
for 7=0.01 and 7=1 show striking variations with respect to |r — s|. For a fixed / (low forces 
with ///c-Cl or ~1) the decay of Prs{t) not only depends on the strength of the contact, but 
also on its location. At f—lpN, Prs{t) for strong and weak bonds display markedly different 
kinetics, with Prs{t) for weak contact (r, s)=(l, 100) and (20,80) decaying to values close 
to zero on 30/is timescale (Figs, la and lb). However, for (r, s)=(40, 60) the asymptotic 
valueof P^,(i^cx))=P;,5«i0.12 (Fig. Ic). Surprisingly, P^,(t)s for (r, s) = (20, 80) and (40,60) 
for strong bonds are qualitatively similar (Fig. lb and Ic) and approach higher values of 
P;,«fti0.42 and 0.85 respectively, whereas P^^^0.2Q for (r, s) = (l, 100) (Fig. la). Thus, the 
rupture of contacts in the middle of the chain is strongly dependent on the coupling between 
the internal modes of the chain and the dynamics of instability due to the applied force. 
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FIG. 2: The Q—Q plots (empty circles) of lifetimes Prs{t) for single contact (r, s)=(l, 100) (x-axis) 
versus lifetimes Pexp{i) (y-axis), for /=40pN and 7=1 (a) and 7=0.01 (b). The rate K in Pexpit) 
is set to the inverse of the average contact lifetime (the maximum likelihood estimate) [7] . 

We refer to this as the "internal motion dominated" (IMD) regime. The stochastic nature 
of the bond rupture kinetics is also reflected in the decay of Prs{t) for strong and weak 
bonds. For (r, s) = (l, 100), Prs(t) decays faster to lower P^^, due to larger chain fluctuations 
at the chain ends compared to contacts (r, s) = (20, 80) and (40,60) in the middle of the 
chain. In contrast to low forces, for ///c>l, Prs{t) for all contacts decay to zero at long 
times, regardless of the contact strength (Figs. Id-lf). However, Prs{t) for strong contact 
(r, s) = (40, 60), which is close to the midpoint of the chain, decays to zero on a much slower 
40/is timescale, compared to 10/is and 1/is timescale for (r, s) = (20, 80) and (1,100). At 
f=40pN, Prs{t) for weak contacts decay to zero on similar time scales (Figs. Id-lf) which 
shows that at large forces rupture kinetics is solely determined by the instability caused 
by the applied tension, with chain dynamics playing relatively minor role. Thus, as f/fc 
increases, the rupture kinetics become increasingly "tension dominated" , i.e. contact rupture 
is due to the applied force. Interestingly, there is a significant plateau in the decay of Prs{t) 
especially for 7=1. The duration of the plateau increases as |r — s| decreases and / increases 
(Figs. Id-lf), and shows that the kinetics of rupture is non-exponential. 

The lifetime distributions, Prs{t) (insets in Fig. 1), show longer tails for the stronger 
contacts; Prs(t) for the interior contact (40, 60) (Figs. Ic and If), does not have a pronounced 
dependence on / compared to Prsit) for contacts (20,80) and (1, 100) (Figs, la, lb, and Id 
and le). In particular, Prs{t) at f=lpN and f=4:0pN agree quantitatively for strong and 
weak (40, 60)-contacts. This implies that the disruption of the interior contacts is mediated 
by internal chain motions, and hence is in the IMD regime even at high /. The influence of 
internal motion of the chain on the contact rupture kinetics is also reflected in the width of 
Prs{t)s which broadens for remote contacts (20,80) and (40,60). 
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Contact lifetime pdf is non- exponential: Forced rupture of noncovalent bonds is often 
described using two-state kinetics, (rs)— >0, which corresponds to the exponential pdf of 
contact hfetimes, Pexp(t)=Kexp [—Kt] where K is rupture rate. The results for Prs(t) in 
Fig. 1 deviate from Pexp{t), especially when ///c<l or when the contact (rs) is in the 
chain interior. The extent of deviation is assessed using the quantile-quantile (Q—Q) plot 
[Zl |8]. A quantile is a number Xp such that 100 xp% of the probability values are <Xp. The 
Q—Q scatterplot compares two sets of probability values: Xp, the quantile for Prs(t), and 
Hp, the quantile for Pexpif)- If the two sets are similarly distributed, the points fall on the 
reference line [7j. The Q—Q plots of the lifetime distribution for the contact (r, s) = (l, 100) 
at f=AOpN and for 7=1 and 7=0.01 show that Prs{t) is larger (smaller) than Pexp{t) in the 
middle range (tail) of the lifetimes (Fig. 2). The non-exponentiality arises because it takes 
time (rtp) for the force-induced tension to propagate to the contact |9]. 

Interplay of timescales: The non-exponential nature of Prs{t) is determined by the in- 
terplay of two timescales, Ttp and the spectrum of relaxation times {rp} {p=l,2, . . . , N). 
These determine the evolution of Ci(t) due to the applied force f, and the broadening of 
the Gaussian density, C2{t), due to the random force F (Eq. (|2])). To reveal the richness 
in the rupture kinetics, we consider Prs{t) for the (r, s)-contacts close to the chain ends 
in the limits f/fc<l, and f/fc>l. For these contacts, the contribution from even Rouse 
modes is negligible |5], and only the slowest mode with relaxation time ri=.^/Ai contributes 



significantly. As a result, Prs{t)=CQ dRrsexp 



, where Cq is the constant. 



-B ,^ 3(R,,-Ci(f))2 
2C2(t)) 

Ci(t)^R^^(l-e~*/^i), C2(t)~(AR^^)2(l-e-2*/n)^ and R^^ and AR^^ denote the equilibrium 
(r, s)-contact vector and fluctuations, respectively. 

1. Weak force, f / fc<l: In this regime, \Ci{t)\<B, and hence the kinetics of Prs(t) is 
dominated by the broadening of C2{t). By evaluating Prs{t) and using the series expansion 
of the error function, we obtain Prs{t)^exp[—B'^/2C2{t)]. At short times, t<^Ti, and we 
can Taylor expand the exponential function in C2{t), which to the first order in t^ri reads 
e^^*/^i^l— 2t/ri. Then, Prs(^)~e~^^^^/* (ci is a constant). At longer times, when t^ri, 
g-2t/n^]^^ and Prs(i)~exp [— 026"^*/"^^] (c2 is a constant). This implies that for a weak 
force, in the long time limit, Prs{t) of the contacts close to the chain ends is described by 
the double-exponential (Gumbel type) distribution of the largest value for the exponential 
density [TU] . 

2. Strong force, f/fc>l: The kinetics of Prs{t) is dominated by the the dynamics of Ci(t). 
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However, the decay of Prsif) is delayed by Ttp. The tension propagation timescale Ttp^Ti. 
a. Tension propagation regime: {t<Ttp). Since |Ci(t)|<i?, the contact hfetime is strongly 
linked to the random motions of the polymer, which results in broadening of the Gaus- 
sian density due to increase in C2(t). Evaluating Prs(t), and assuming |Ci(t)|~-B<C2(t), 
we obtain Pj.s(t)~(t/ri)^e~'^^(*/'^^^^''^ (03 is a constant). Thus, Prs(^) scales with force as 
-frs(^)~/^e~^'*-^ (c4 is a constant), b. Kinetic regime: (t>Ttp). The contact lifetime is still 
linked to the random motions of the polymer (broadening of C2{t)) but is dominated by the 
dynamics of Ci{t)>B. In this regime, e"*/'^!^! in Ci(t) and e"^*/^i<l in C2(t), and Prs(^) 
is again double-exponential (Gumbel) density [ID], i.e. Prs(^)~exp [—056^^*/^^], where C5 is 
a constant. We find Prs{t)^pe~'^^'^ , where cq is a constant. 

The contact lifetimes for two contacts: The Hamiltonian of the chain with two con- 

tacts (r, s) and (r', s') is given by H='^Y.n=2 R^,„-i+fRm+ Y.i^,)W^^B,,-\^i^)%^ 
where hij {Bij) are the equilibrium (critical rupture) distance, which we rewrite as 
H=Hrs+Hr's', where Hij='^^-^YlK,n-iH^Ni+^i^Q{Bij-\Rij\)R]y The population 
is P(t)=Prsys'(t)+Poys'(t)+Prsfl(t): whcrc Prs,r's'(t)=Prs(t)Pr's'(t) dcuotcs that both con- 
tacts are intact, Poys'(t)=Prs~>o{t)Pr's'(t) means that (r, s) is ruptured but {r',s') is in- 
tact, and Prsfiit)=Prsit)Pr's'->Qit) implies that {r',s') is disrupted but (r, s) is intact. 
Pij^o(t)=^~Pij(t) is the population of conformations with the disrupted contact {i,j). 
Using r'^r, and brs^br's'=b, the Hamiltonian for two nearest neighbor contacts 

is H=^^Yl^l,n-i+^^Ni +^e{B-\Rrs\)Rl„ whcrc b=b/V2 is the rescaled distance. 
Ptotit) iptotit)) is computed using Prsif) {prsif)) for single contact (r, s). 

Plots of Ptot{t) for two separated strong contacts at f=40pN and 80pA^ show that, in 
general, increasing / facilitates more rapid rupture of the binary contacts (Figs. 3a-3c). 
The presence of interior contacts delays their rupture (compare Figs. 3a with 3b and 3c) 
even though ///c>l. The time dependence of Ptotif) for two separated contacts (Figs 3a-c), 
after a lag phase, can be analyzed using two distinct exponential functions, which shows 
the rupture occurs by sequential kinetics, with {(r, s), (r', s']}^{(0), (r', s')}^{(0), (0)]}. 
The rate limiting second step is the rupture of the more interior contact. Comparison of 
Ptotit), for 7=1 and f=AOpN, for the nearest neighbor contact pairs (r, s) = (l, 100) and 
(r',s') = (2,99), (20,80) and (21,79), and (40,60) and (41,59) (Fig. 3d-3f) with P^,(t) for 
single contact (r,s) = (l, 100), (20,80), and (40,60) (Fig. Id-lf) shows that the double bond 
increases stability of the contact, especially for the interior contact (40,60). The chain 
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FIG. 3: The survival probabilities Ptot{t) for contact pairs (1,100) and (20,80) (a), (1,100) and 
(40,60) (b), and (20,80) and (40,60) (c), (1,100) and (2,99) (d), (20,80) and (21,79) (e), and 
(40,60) and (41,59) (f), for /=40pN (solid lines) and 80pN (dashed lines), all with 7=1. The life- 
time distributions ptot{t) are shown in the insets. Time scales vary greatly in all panels. Structures 
are for the intact pair (20,80) (40,60), disrupted contact (20,80) and intact contact (40,60), and 
both contacts disrupted are shown. 

with double bonds ruptures in a single step, {(r, s), (r', s')}^{(0), (0)}, i.e. the rip occurs 
cooperatively, so that Ptot{t) decays with a single rate constant at long times. The decay 
of Ptot{t) for separated contacts, (r, ,s) = (20, 80) and (r', s')=:(40, 60), is slower than Ptot{t) 
for the nearest neighbor contacts, (1, 100) and (2,99) (Figs. 3c and 3d), and P/^^ is larger, 
which shows that the persistence of binary interactions critically depends on the location 
of contacts in the chain. The lifetime distributions show that increasing the force from 
f—AQpN to f=80pN results in the decreased stability of contacts and shorter lifetimes (Fig. 
3). The contact lifetimes increase when they are located in the chain interior (Fig. 3). The 
lifetimes of the nearest neighbor contacts, (1, 100) and (2,99) (Figs. 3d and 3e), are shorter 
than that for the pair, (40,60) and (41,59) (Fig. 3f), which shows the importance of the 
location of binary contacts with respect to the point of application of the force. 

Conclusions: The present work shows that the shape of the free energy landscape can 
only be discerned by analyzing the topology-dependent lifetime distributions of adhesive 
contacts that reflect the structure the molecules. The continued increase in the temporal 
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resolution in laser optical tweezer experiments should allow the prediction that, the internal 
motions of biomolecules are coupled to global fluctuations, to be tested. The present theory 
provides a conceptual framework for interpreting such experiments [T]. In contrast, the 
rich behavior predicted here cannot be obtained by solving for the dynamics in the exactly 
calculable equilibrium free energy profile in the presence of non-zero force. Accounting for 
the non-exponential distribution of unbinding lifetimes will require models that reflect the 
interplay between local chain motions and tension-induced global unfolding. 
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